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The observation of a novel type of commensurability resonance in two-dimensional, hexagonal 
antidot lattices is reported. These resonances have a classical character and occur at magnetic 
fields above the resonance that corresponds to the cyclotron motion around a single antidot. The 
resonances are visible only for antidots with effective diameters larger than 50 % of the lattice 
constant. Simulations reveal that they originate from quasi-stable electron trajectories that bounce 
between three neighboring antidots. This interpretation is backed by the observation of large-period 
Aharonov-Bohm type oscillations at low temperatures. 



PACS numbers: 73.23.-b, 73.63.-b 



INTRODUCTION 



The study of electron transport in artificial, two- 
dimensional periodic potentials has revealed a variety 
of interesting phenomena over the past 15 years. One 
variant of such systems are antidot lattices, i.e. periodic 
potentials with maxima above the Fermi energy of the 
two-dimensional electron gas (2DEG) 0, 0, Most 
strikingly resonances in the longitudinal magnetore- 
sistance are found, which can be interpreted in terms 
of classical cyclotron orbits that are commensurate 
with the antidot lattice Q- A more thorough classical 
treatment based on the Kubo formalism Q| has revealed 
that these resonances actually have their origin in 
the magnetic field dependent mixture of chaotic and 
regular trajectories, where the latter remain pinned in 
weak electric fields. Moreover, this theory provides an 
explanation for the observation of a negative Hall effect 
in weak magnetic fields B 5]. Another striking effect 
is the occurrence of _B-periodic oscillations, which have 
been explained within a semiclassical theory by a few yet 
dominant quantized periodic orbits |g, |2| ■ More recently, 
signatures of the famous fractal energy spectrum of such 
potentials, also known as the Hofstadter butterfly, have 
been observed 

A large fraction of these studies was performed on square 
[1 Sll II H H E3 or rectangular lattices, while 
only a few ex peri ments on hexagonal lattices have been 
published [Tl 111 HI [TEj . For the majority of the effects, 
the lattice type is, in principle, irrelevant. Nevertheless, 
hexagonal lattices show some peculiarities that are 
absent in other types of antidot lattices. In particular, 
Altshulcr-Aronov-Spivak oscillations have been observed 
around B = in hexagonal lattices 0, while 
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Aharonov-Bohm oscillations can be detected at larger 
magnetic fields 0, 0. Also, it has been suggested 
recently that scattering centers with a short range 
hexagonal order may be responsible for a phenomenol- 
ogy resembling that one observed in metal-insulator 
transitions in two dimensions 0|. Moreover, a detailed 
understanding of the transport in two-dimensional 
hexagonal lattices is of broad relevance, since this type 
of lattice forms via self-organization on a mesoscopic 
scale in a variety of systems, such as vortex lattices in 
type II - superconductors 01 or diblock copolymers • 



Here, we report the observation of a novel type 
of magnetotransport resonance in hexagonal antidot 
lattices with large effective antidot diameters d e > 0.5a, 
where a denotes the lattice constant. These resonances 
occur at cyclotron radii smaller than the antidot radius. 
Their signatures are visible only in the longitudinal 
component of the magneto-resistivity tensor, while tem- 
perature dependent measurements indicate a classical 
origin. Furthermore, Aharonov-Bohm type oscillations 
are observed at low temperatures that correspond to 
enclosed areas much smaller than the size of a lattice 
unit cell. These resonances are therefore attributed to 
quasi-stable trajectories which are localized in between 
three neighboring antidots. Numerical simulations of 
the magnetoresistivity based on the Kubo formalism 
support this interpretation and allow us to specify the 
characteristic trajectories. 



The paper is organized as follows. In Section II, the 
sample preparation and the experimental setup is dis- 
cussed. Section III is devoted to the experimental char- 
acterization of the resonances. In Section IV, we present 
numerical simulations, identify the relevant trajectories 
and relate these results to the experiments. A summary 
and conclusion is given in Section V. 
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II. SAMPLE PREPARATION AND 
EXPERIMENTAL SETUP 

Conventional modulation doped GaAs / AIq^Gclq^As 
- heterostructures with a two-dimensional electron gas 
(2DEG) 83nm below the surface have been grown by 
molecular beam epitaxy. The 2DEG has an electron 
density of n = 2.5 x 10 15 m -2 and a carrier mobility of 
90m 2 /Vs, corresponding to an elastic mean free path of 
7.6/im and to a Drude scattering time of r = 35ps at liq- 
uid helium temperatures. A standard Hall bar geometry 
has been defined by optical lithography and wet chemical 
etching. The 2DEG is accessed electrically via Ni/AuGe 
Ohmic contacts. The hexagonal antidot arrays have been 
patterned by electron beam lithography and subsequent 
ion beam etching with low-energy (350eV)Ar + ions [l9| . 
The etch depth was 50nm. All antidot arrays have a lat- 
tice constant of 600nm, while lithographic antidot diame- 
ters of di = lOOnm, 200nm, 250nm and 320nm have been 
patterned. The measured array area was 50/im x 25/im. 

Transport experiments have been carried out in a he- 
lium gas flow cryostat with a variable temperature insert 
and a superconducting magnet with a maximum field of 
B = 8T, which is applied perpendicular to the plane of 
the 2DEG. The temperature has been varied from 1.6K 
to 24K with an accuracy of 0.1K. For measurements 
at lower temperatures, a top loading dilution refrigera- 
tor with a base temperature of 80mK has been used. An 
AC current (frequency 33Hz, amplitude 5nA) was passed 
through the antidot array. The longitudinal (p xx ) and 
Hall (p X y) components of the resistivity tensor have been 
determined from voltage differences measured in the cor- 
responding four- probe setups. 



III. EXPERIMENTAL RESULTS 

Fig. 1 shows the measurements of 

Pxx & 

function 

of B for arrays with lithographic antidot diameters of 
di = lOOnm, 200nm and 320nm. Due to lateral depletion 
of the 2DEG around the antidots of about lOOnm, their 
electronic diameters d e are larger. They can be estimated 
from the Shubnikov-de Haas oscillations observed in the 
antidot arrays, which give the Fermi wavelength in the 
array as well as, via their onset, the number of occupied 
transverse modes at the bottleneck formed by adjacent 
antidots. This way, effective electronic diameters of d e « 
300nm, 400nm and 550nm, respectively, are found. 

In the array with the smallest antidots (Fig. 1 (a)), 
typical commensurability oscillations in p xx are observed 
that, within the simple commensurability picture, in 
which the cyclotron diameter r c = m*VF / (eJ3) matches 
the lattice constant in the resistivity maxima, can be at- 
tributed to orbits around 1, 3, and 7 antidots [lj. Here, 
the resonance around one antidot occurs at B = 270mT. 
In comparison to the decay of the Shubnikov - de Haas os- 
cillations with increasing temperature, the temperature 
dependence of the commensurability oscillations is weak, 
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FIG. 1: Longitudinal magnetoresistivity of the hexagonal 
antidot arrays at temperatures of T — 1.5K, 2.5K, 4K, 8K, 
12K, 16K, and 24K, as observed for arrays with a = 600nm 
and different values of di (inset). In (a), the well known res- 
onances corresponding to cyclotron orbits around 1, 3, and 7 
antidots are observed. In (b) and (c), a known resonance is 
observed around B = 130mT. The additional resonances are 
indicated by arrows. The inset in (c) shows a zoom-in from 
the main figure of the region of interest. 
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FIG. 2: Hall resistivities of the 3 arrays of Fig. 1, measured 
at T = 1.5K. In contrast to the conventional commensura- 
bility oscillations, the additional resonances cannot be identi- 
fied. Note the negative Hall effect in the array with the small 
antidots (inset). 
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FIG. 3: Magnetoresistivity of the array with di — 320nm 
as measured in the dilution refrigerator in the magnetic field 
regime of the novel resonances. A smooth background has 
been subtracted. Here, the quasi-periodic Aharonov-Bohm 
oscillations are modulated with a period of SB m I20mT. 



which indicates their classical origin. Note that there 
is an additional shoulder at B = 200mT and a peak at 
70mT with a weak temperature dependence (vertical ar- 
rows in Fig. 1(a)), which cannot be identified in this 
commensurability picture. 

Fig. 1 (b) shows the corresponding measurements on 
a sample with di = 200nm. Here, only one strong reso- 
nance with a classical character is visible at B = 130mT, 
which from simulations as described below can be at- 
tributed to a superposition of the commensurability os- 
cillations identified in Fig. 1 (a). In addition, a weak 
resonance is observed at B = 400mT. It is clearly not 
related to the Shubnikov - de Haas oscillations that set in 
around B = 0.6T, and shows a weak temperature depen- 
dence. Such a structure is also present in an array with 
di = 250nm (not shown). Moreover, in Fig. 1 (b), a small 
feature is observed between B = 0.5T and 0.7T which, 
however, cannot be clearly distinguished from Shubnikov 
- de Haas oscillations. In arrays with even larger an- 
tidot diameters (Fig. 1 (c)), three maxima, superim- 
posed on a Shubnikov - de Haas resonance, are detected 
at B « 0.45T, 0.87T, and 1.15T. The main focus of the 
present paper is on the investigation of these resonances. 
The corresponding Hall measurements are shown in Fig. 
2. While the structure of the conventional commensura- 
bility oscillations is reflected in p xy (B), no signatures of 
the novel resonances are detected here. 

To obtain further insight into the character of these 
structures, we have studied the magnetotransport of the 
array with di = 320nm in a dilution refrigerator, see 
Fig. 3. Around B = 0, the well-known Altshuler- 
Aronov-Spivak oscillations with a period of SB = 8mT 
are observed (not shown) |13|. which evolve into quasi- 
periodic Aharonov-Bohm type oscillations with a period 
of 5B = 16mT for B > 60mT. The Aharonov-Bohm pe- 
riod corresponds to a characteristic area of A = h/eSB = 



2.6 x 10~ 13 m 2 , in rough agreement with the area of a lat- 
tice unit cell (A ce u = 3.1 x 10~ 13 m 2 ). This evolution has 
been reported and interpreted previously 0]. An addi- 
tional, superimposed oscillation, however, with a period 
of SB f» 120mT is observed in the magnetic field regime 
where the novel classical resonances are detected. Inter- 
preting this oscillation in terms of an Aharonov-Bohm 
effect, the enclosed area equals A « 3.4 x 10~ 14 to 2 , which 
corresponds to about 11% of a unit cell. All periodic os- 
cillations smear out as the temperature is increased, with 
the amplitude decaying roughly oc 1/T; the oscillation 
with the larger period persists to higher temperatures. 



IV. INTERPRETATION AND DISCUSSION 

We interpret our results in terms of commensurabil- 
ity oscillations that exist in the open electron pockets 
formed in between three adjacent antidots. Quasi-stable 
orbits may form by consecutive reflections of the elec- 
trons at the antidot walls. It seems plausible that such 
trajectories become more stable as d e /a increases. We 
determine the characteristic area of such a pocket from 
Apocket = \{A cei i - Ti{d e /2) 2 ) and estimate in our large 
antidot array to A poc k e t ~ 3.6 x 10 _14 to 2 , i.e., signifi- 
cantly smaller than the area of the unit cell. Hence, an 
increased Aharonov-Bohm period can be expected, al- 
though the exact period depends on the details of the 
closed trajectory. 

In order to substantiate this interpretation, we have 
performed model calculations based on the Kubo for- 
malism J2jj. The antidot lattice is modelled by a two- 
dimensional, hexagonal array of hard- wall cylinders of 
the estimated electronic diameter d e for the samples of 
Fig. 1. For each value of B, 10 5 electrons with a fixed 
Fermi energy adapted to the experimental conditions 
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(Ep = 8.86meV) are injected at random positions and 
with random velocity directions within a unit cell of the 
lattice. Their trajectories are calculated to a length of 
50/xto, corresponding to a time of flight of 240ps. The 
electron velocity correlation function (vi(t, B)vj(0)) is 
calculated, where the brackets denote averaging over all 
trajectories [2l|. Within the Kubo formalism, the com- 
ponents of the magneto-conductivity tensor for a degen- 
erate two-dimensional electron gas are obtained from 

oo 

°«( B ) = 5"/ {ui{t,B)vM)z- tlT at (i) 



In Eq. (JJJ, to* = 0.067TO e is the effective electron mass 
in GaAs. Elastic scattering due to random impurities is 
taken into account via the exponential cutoff function 
that appears in the integrand. Here, a Drude scattering 
time of t = 35ps has been chosen, in accordance with the 
mobility of the pristine 2DEG. This model is a simplified 
version of the more thorough treatment presented in Ref. 
0. In particular, neither the driving electric field nor the 
finite slope of the antidot walls are taken into account. 
Nevertheless, our simplified model allows us to identify 
the trajectories that generate the additional resonances. 
In Fig. 4, the results of the simulated p xx (B) for the 
three arrays under study are reproduced. For the largest 
antidots (d e — 550nm), pronounced resistivity minima 
are found around B = 0.78T, 1.4T, and 1.85T, while the 
commensurability resonance around a single antidot is 
barely visible at B ss 0.3T. 

To gain more insight into the origin of these structures, 
we have taken Poincare sections along the perpendic- 
ular bisector of the line connecting the centers of two 
neighboring antidots, see Fig. 5. As the magnetic field 
increases to B = 0.6T, stable regions emerge with a 
structure as shown in Fig. 5 (a). Around the resistivity 
maximum at IT, this pattern fades and evolves into a 
second one shown in Fig. 5(b), which is most pronounced 
around the minimum in p xx at 1.4T. As one goes to 
larger magnetic fields, the stable regions gradually fade 
until they can no longer be identified. We pick the 
Poincare sections at B = 0.6T and 1.2T for a detailed 
discussion. In both cases, stable regions are found which 
contain the two types of trajectories shown in the right 
parts of Figs. 5 (a) and (b), respectively. The first type 
is labelled by 1 and consists of skipping orbits around an 
individual antidot. This type is present for all magnetic 
fields above 0.5T and covers a continuously increasing 
area in the Poincare section as B is increased up to 
2T. Such skipping orbit trajectories have already been 
discussed for square antidot lattices ^(J- The second 
type are triangularly or rosette-type shaped trajectories 
that form quasi-closed figures by consecutive reflections 
at walls of adjacent antidots. They are labelled by 2 - 4 
in Fig. 5. The areas of the corresponding regular islands 
in the Poincare sections have broad peaks, centered 
around the minima in p xx at B = 0.78T and 1.4T. 
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FIG. 4: Simulated longitudinal magnetoresistivity of the 
array with d e = 550nm. Weak resonances are found in 
0.6T < B < 2.2T. The inset shows the corresponding simula- 
tions for the arrays with d e = 300nm and d e = 400nm. Here, 
the star denotes an unidentified resonances, while the cartoon 
to the upper right sketches the type of trajectory found in the 
feature indicated by the arrow. 



Type 2 - trajectories with two consecutive bounces at 
the same wall are present only for B > IT. Hence, we 
attribute the structure in the resistivity of the array 
with the large antidots to the triangular or rosette - 
shaped trajectories. They are most important around 
the minima of p xx , which is a consequence of the large 
magnetic field, since localized states that generate a 
minimum in a xx {B) translate into a minimum in p xx {B) 
for p xy 3> p xx , as is the case here. Thus, a characteristic 
feature of the novel commensurability resonances which 
distinguishes them from those known previously is the 
fact that the magnetic field at which the quasi-stable 
trajectories are most pronounced corresponds to a 
minimum in p xx . 

Despite this qualitative agreement between the exper- 
imentally observed features and the simulations, quan- 
titative differences are apparent. In particular, the po- 
sitions of the measured minima and maxima in p xx dif- 
fer from the simulated ones, and their amplitudes are 
weaker, an effect that, due to the observed temperature 
independence, cannot be attributed to thermal smearing. 
Also, the structure stemming from rosette-type trajecto- 
ries can be identified in simulations of the Hall resistiv- 
ity (not shown), in contrast to our Hall measurements. 
One possible reason for these discrepancies may be non- 
specular scattering at the antidot walls. Low-energy ion 
beam milling has been the technique of choice for pat- 
terning of the antidots, since it is established as a method 
to generate nanostructures with walls of high specular- 
ity around 85% [l^. We expect that the small fraction 
of non-specular scattering may lead to some degree of 
smearing, but most likely not to a displacement of the 
extremal points. Rather, we interpret these deviations 
as a manifestation of soft antidot walls, probably in com- 
bination with some disorder in the antidot sizes, shapes 
and positions, as well as the driving electric field. We 
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are not aware of experimental work that determines the 
wall steepness for the etching technique we used. How- 
ever, it seems plausible to assume that the steepness is 
comparable to that one obtained in Ga[Al]As samples of 
similar electron density, patterned by other dry etching 
techniques, like those used in Ref. [2|. More elaborate 
simulations that use the wall steepness as a parameter 
0, have been able to establish a close to perfect agree- 
ment with the experimental observations of Ref. . As 
an example, we note that a peak is observed in the exper- 
iments on antidot arrays with small diameters at small 
magnetic fields (Fig. 1 (a)), which cannot be attributed 
to a trajectory by our simulations. A similar discrepancy 
in square lattices 0] could be resolved by the softness of 
the walls a fact which suggests that the deviations in 
our study can be explained by soft walls as well. In gen- 
eral, soft walls tend to deform the trajectories, which in 
turn causes some resonances to be significantly displaced, 
while the sensitivity to the wall steepness depends on the 
resonance Q . Such a deformation would be in tune with 
the observation that the areas enclosed by the simulated 
trajectories (Fig. 5) are smaller than the measured en- 
closed areas (Fig. 3). Also, the peak around 0.35T in 
the simulated array with d e = 300nm can be identified 
in a Poincare section of the hard wall simulation (Inset 
in Fig. 4), but is not observed experimentally. It seems 
plausible to us that such a trajectory is very sensitive 
to an even moderate softening of the wall. To study 
whether more realistic simulations, or experiments on 
samples with steeper walls [lQj, are able to resolve the 
discrepancies, is beyond the present study. 
For the sample with the smallest antidot diameter, the 
simulations reproduce the experiments remarkably well 
in the range r c > a (Figs 1 (a), 2 and 4). The mea- 
sured negative Hall (Fig. 2) effect is reproduced in our 
hard wall simulations (not shown), in contrast to hard 
wall simulations performed for square lattices 0]. In 
square lattices consisting of antidots with similar ratios 
d/a, simulations do not suggest the existence of type 2 - 
orbits, independent of the wall steepness Thus, the 
wall steepness seems less critical in hexagonal arrays, but 
becomes more relevant as d e /a increases. We therefore 
conclude that a hexagonal lattice strongly favors the for- 
mation of stable orbits in between antidots as compared 
to square lattices. 



V. SUMMARY AND CONCLUSION 

In summary, we observe novel resonances in hexagonal 
antidot lattices with large antidot diameters. Classical 
simulations suggest that they can be attributed to 
classical, quasi-periodic trajectories of various shapes 
that are caught in between three adjacent antidots. 
This interpretation is supported by measurements 
at temperatures below IK, where the samples show 
Aharonov-Bohm type oscillations with a characteristic 
area similar to that one of an electron pocket in between 



(a) B = 0.6T 




-1 1 
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FIG. 5: Poincare - sections (left) taken from the array with 
d e = 550nm at B = 0.6T (a) and 1.2T (b). The sections are 
taken along the y-direction at x = 0, see the right part of (a), 
where the gray areas denote the antidots. The characteristic 
trajectories (right) found in the stable regions are assigned 
according to the numbers in the stable regions, see text. Note 
that only the intersections of the trajectories with v x > are 
plotted in the Poincare - sections. 

antidots. It remains to be seen whether such resonances 
are also visible in other types of antidot lattices with 
optimized parameters. In order to improve the quantita- 
tive agreement with simulations, it would be interesting 
to perform more sophisticated calculations with realistic 
potential landscapes, inclusion of the driving electric 
field, disorder, and the wave nature of the electrons. 
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